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ANALYSIS OF THE DYNAMIC RESPONSE OF A SUPERSONIC INLET 
TO FLOW-FIELD PERTURBATIONS UPSTREAM 
OF THE NORMAL SHOCK 
by Gary L. Cole and Ross G. Willoh 
Lewis Research Center 

SUMMARY 

A linearized mathematical analysis of supersonic inlet dynamics is presented. 
Attention is concentrated on determining the response of normal shock position and sub- 
sonic duct pressures to flow-field perturbations upstream of the normal shock in mixed- 
compression inlets. The analysis is based on a previous NASA report which dealt pri- 
marily with perturbations downstream of the normal shock. 

The inlet duct cross-sectional area variation is approximated by constant area sec- 
tions. This approximation, in combination with a linearized analysis, results in one- 
dimensional wave equations for each duct section. The supersonic and subsonic flow re- 
gions are separated by a movable normal shock. A choked exit is assumed for the inlet 
exit condition. The analysis leads to a closed-form matrix solution for the shock posi- 
tion and pressure transfer functions. 

The analysis was compared on a frequency response basis with a method-of- 
characteristics solution. The agreement in both amplitude ratio and phase was excellent. 

Analytical frequency response results were also compared with experimental data. 
The phase angle results were generally in good agreement. Amplitude ratio response 
curves, although not in as good agreement as phase data, showed generally good agree- 
ment in shape. Some shifts in low-frequency gain were found, however. 

INTRODUCTION 

Propulsion system performance of supersonic aircraft depends greatly on the type of 
inlet system being used. Mixed- compression inlets are best for achieving high propul- 
sion system performance in aircraft that fly at Mach numbers in excess of about 2.0. 
Because of the nature of mixed- compression inlet design, a normal shock wave usually 
exists within the inlet during supersonic operation. For best inlet performance the nor- 
mal shock must be positioned near the throat where the shock is susceptible to being 



displaced by disturbances arising from such things as atmospheric perturbations, air- 
craft maneuvers, and changes in engine operation. A downstream displacement of the 
shock results in a loss of inlet performance. K the shock is displaced in the upstream 
direction, it is in danger of being expelled from the inlet. This is referred to as inlet 
unstart. Inlet unstart may result in undesirable consequences such as compressor stall, 
combustor flameout, reduced propulsion system thrust and increased vehicle drag. To 
counteract such possibilities, mixed-compression inlets are provided with variable 
geometry features that can be automatically controlled to keep the shock at the desired 
position. The design of these shock-position control systems requires knowledge of the 
shock-position dynamic response to perturbations in the inlet. The purpose of this 
report is to develop an approximate mathematical analysis for predicting the dynamic 
response of shock position and pressures, in mixed-compression inlets to perturbations 
upstream of the normal shock. 

In the past a great deal of attention has been given to controlling shock position 
against airflow perturbations originating downstream of the normal shock (refs. 1 to 6). 
Normal shock motion due to perturbations downstream of the shock has been examined 
analytically (e.g., refs. 7 and 8). Reference 9 presents a mathematical analysis that is 
valid for obtaining shock responses to perturbations downstream or just upstream of the 
normal shock. However, the upstream terms were dropped and consideration was 
restricted to the downstream case. Reference 10 does deal with upstream perturbations 
and presents an analysis that includes both storage-volume and Helmholtz-mass effects. 
Reference 11 presents transfer functions for shock position to upstream perturbations. 
The transfer functions vary in complexity from first to fourth order, and results are 
compared with a method-of-characteristics solution. The fourth-order model is a lin- 
earized version of the analysis in reference 10. 

In this report a mathematical analysis is presented that combines a set of linearized 
equations across the normal shock with an exact solution of the linearized one- 
dimensional wave equation. This analysis is based on the analysis of reference 9. Use 
of this technique avoids the complexity of the method of characteristics while still pre- 
dicting the resonances of a distributed parameter system. However, flow -field discon- 
tinuities due to oblique shock waves in the supersonic duct are neglected. The analysis 
presented is more exact than conventional lumped-parameter techniques and for fre- 
quency response results no more complicated in application. If transient responses are 
desired the equations are especially suitable for analog simulation. The analysis is fol- 
lowed by a matrix solution that provides a simple means for obtaining frequency 
responses on a digital computer. Finally, analytical frequency responses are compared 
with results from a method-of-characteristics solution found in reference 11 and with 
experimental data obtained during a wind tunnel program. Dr. Frank Barry of the 
Hamilton Standard Company supplied us with information regarding the method-of- 
characteristics analysis in addition to that in reference 11. 
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MATHEMATICAL ANALYSIS 


The analysis presented in this report is an extension of the analysis in reference 9. 
In this report prime consideration is given to perturbations upstream of the normal 
shock. The symbols used in this report (and defined in the appendix) are consistent with 
those of reference 9. 

A schematic of an ideal mixed-compression inlet is shown in figure 1. In refer- 
ence 9 a set of linearized equations relating normal shock position to adjacent param- 
eters was developed using total pressure, flow, and entropy as state variables. In addi- 
tion, a linearization of the compressible flow equations written in terms of the state 
variables was used to develop a set of wave equations for the subsonic duct, hi the anal- 
ysis of this report the same wave equations are applied to both the subsonic and super- 
sonic flow regions. The normal shock position equations are used to couple the subsonic 
and supersonic duct sections and a choked exit is used as the downstream boundary con- 
dition for the subsonic duct. 


Shock- Position Dynamics 

Equations relating normal shock position and velocity to adjacent upstream and 
downstream variables (total pressure, airflow, and entropy) were obtained using the fol- 
lowing assumptions: perfect-gas relations, a constant specific -heat ratio of 1.4, and 
negligible change in duct area across the normal shock. The normal shock equations are 
(ref. 9, eqs. (1) to (3)): 
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where the subscripts 1 and 2 designate fixed stations just upstream and just downstream 
of the shock (fig. 1), respectively. There are two differences between equations (l)to (3) 
herein and those of reference 9. One is that shock position is normalized by the radius 
of the inlet cowl lip R c instead of by the length of the subsonic duct (Lrp in ref. 9). (R c 
is a more convenient normalizing parameter because it is a constant, whereas L^, de- 
pends on the shock operating point. ) The other difference is that equations (1) to (3) herein 
have the quantity Wj; in reference 9 the subscript 1 was omitted. In this analysis the 


airflow rates will be normalized by the steady-state airflow at the most upstream sta- 
tion, designated as W M . To obtain the correct normalized airflow rates the AW^/Wj 
terms in equations (1) and (3) must be multiplied by the ratio W^/W^, and equation (2) 
must be multiplied by Wj/W m . 

Equations (1) and (3) contain the variable Aa T 1 /a T , which is not independent of 
(ASj/R) and (AP t j/FVp j)- The Aa T j/a T terms can be eliminated by means of the 
following equation (ref. 11): 


Aa T _ 1 / AP T + AS 
&t 7 [p T R 


(4) 


Using equation (4) to eliminate Aa-p terms and including the multiplications 

involving W< and W^, equations (1) to (3) become 
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Equations (5) to (7) can be written in terms of the Laplace transform as 
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Supersonic inlets often have provisions for bleed flow in the vicinity of the normal 
shock operating point to increase inlet stability. The analysis of reference 9 did not 
include the effects of stability bleed in the normal shock equations (8) to (10). If the 
effects of small perturbations in pressure are neglected, such bleed flows are propor- 
tional to shock position. Reference 12 shows that for such a bleed equation (9) becomes 

*2 = *1 + ( C B + c r c 3 s ) *s (ID 

where 

(aV^) 

B (AX a /R c ) 

Duct Wave Equations 

The wave equations were derived in reference 9 as the solution for the subsonic 
duct. The following assumptions were made: perfect gas relations, a constant specific - 
heat ratio of 1.4, an inviscid nonheat-conducting fluid, and negligible change in duct area 
that might result from boundary-layer effects or changes in inlet geometry. 

These assumptions are not restrictive to subsonic flow. However, discontinuities 
in entropy and total pressure due to oblique shock waves in the supersonic portion of the 
duct are neglected. The linearized one -dimensional wave equations are (eqs. (96) 
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to (98), ref. 9): 
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The solutions of equations (12) to (14) may be used to relate properties between two 
stations, n and n + 1, in a duct of length L. In Laplace transform notation the 
results are 
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Equations (15) and (17) represent the propagation of waves traveling downstream. 
Equation (16) also represents the propagation of a wave, but the direction depends on the 
value of M. In supersonic flow (M > 1) the wave travels downstream and in subsonic 
flow (M < 1) it travels upstream. 

Equations (15) to (17) are valid only for a constant-area duct. Variations in inlet 
duct area are approximated by dividing the inlet into a series of constant-area ducts. 

The Mach number in a duct section is averaged by summing one-half of the values at 
each end of the duct section. The computation of a, y and a, r, 6 for each 

section is based on the average Mach number and speed of sound in the section. 

Equations (8), (10), (11) and (15) to (17) account for shock position and duct dynam- 
ics. The set provides a basic mathematical model that can be applied to a mixed- 
compression inlet. For the solution of a specific problem, bleeds and bypasses must be 
accounted for, an equation for the termination of the final subsonic section must be pro- 
vided and the disturbance must be related to the system variables (P, and S). 

Disturbance Equations 

When perturbations occur ahead of the normal shock (in the supersonic flow region), 
the perturbation actually represents the upstream boundary condition. This is true 
because the waves represented by equations (15) to (17) all travel in the downstream 
direction. Thus for a perturbation 6 at station 0 upstream of the inlet 

p 0 = qi & 

W 0 = q 2 6 ) (18) 

S 0 = q 3 6 ^ 

where the q's represent constant coefficients. A pure flow perturbation, for example, 
would have q^ = 0, q 2 = 1, and q 3 = 0. 
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Exit Boundary Condition 


As in reference 9 it is assumed that the airflow at the inlet exit is choked. With the 
additional assumptions that the specific heat ratio is constant and equal to 1.4, the equa- 
tion for the exit boundary condition is (ref. 9): 
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If the inlet has bleed or bypass flows, the steady- state exit flow W E will not be 
equal to the inlet flow W M . Since all flows are normalized to W^, equation (19) must 
be multiplied by the ratio W^W^. Therefore, 
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Assuming that the exit area is constant (AAg = 0) and taking the Laplace transform, 
equation (20) becomes 
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W. 
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( 21 ) 


It should be noted that any boundary condition that can be described in terms of the sys- 
tem variables (P, W, S) could be used. 


Bleed and Bypass Airflow Equation 

Supersonic inlets often have bleed and bypass airflows for such things as boundary- 
layer control and matching inlet airflow to engine airflow requirements. These airflows 
(hereinafter called bypass flow) should be accounted for because they affect the gains of 
normal shock position and inlet pressures to the perturbation variable. Equations for 
these airflows will now be derived based on choked flow conditions in order to be consist- 
ent with the experimental inlet that will be discussed later. A typical bypass flow in the 
wall of an inlet duct is shown in figure 2. It is assumed that the bypass flow occurs over 
a zero length section of duct and that the bypass area is constant. 

Since the bypass is choked, the form of equation (21) applies such that 
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Since the bypass is in the side of the duct, it is assumed that the total pressure 
acting on the bypass is equal to the static pressure at station n. Hence, P^y = p n 
where p = (AP g (s)/P g ) and in addition §^y = ^ n - Reference 9 showed that 
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Therefore, in terms of variables at station n, equation (22) becomes 
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MATRIX SOLUTION OF EQUATIONS 


Even though time domain solutions of the complete set of equations is difficult be- 
cause of the algebraic complexity, a closed-form matrix solution for the shock posi- 
tion transfer function can be obtained. In matrix form the normal shock equations (8), 
(10) and (11) become 
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or to use a short notation form 


T 2 = HT 1 + GX S 


(25) 


In general, to relate conditions at two stations n and n + 1, the wave equations 
(15) to (17) are used. In matrix form they become 
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(where B“* is the inverse matrix of B). 

It was shown in reference 9 that the elements of E, e- are 
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The equations relating variables across a choked bleed or bypass are 


P n + 1 = P n 


W n + 1 = W n " W by 


S , = S 
n+1 n 


Using equation (24) for W^y > the matrix form of these equations becomes 
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Equation (21) for the exit boundary condition in matrix form is 
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All that remains is to write a matrix equation relating the perturbation and the inlet 
variables. In matrix form the expressions describing the perturbation (upstream bound- 
ary condition), equations (18) become 
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33 


or 

T q = Q6 (31) 

The equations ((25), (27), and (29) to (31)) can then be used to form a matrix solution 
for determining the transfer function between shock position and the perturbation. To 
demonstrate the matrix solution, an example will now be presented. 

Consider a case in which a perturbation in total pressure occurs at the inlet spike- 
tip station. The inlet for the example is shown schematically in figure 3. To simplify 
the example, only four sections were used, but more sections could be used to better 
approximate the variation in duct area. The analysis begins with equation (25) 

T 2 = HTj + GX g (25) 

Variables at station 1 can be related to those at the cowl-lip station, O', by using 
equation (27) 
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Substituting for Tp equation (25) becomes 


T 2 = HE b E a T 0 + GX s 


Equation (31) relates Tq to the disturbance and is repeated here: 

T 0 = Q5 

where, for this example, 


Q = 


and 6 = P r 


since P-j«,0 is the only variable perturbed. Using (31), equation (32) becomes 

T 2 ' HE bW> + G * S 

Using equations (27), (29), and (30) results in 
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Therefore, multiplying both sides of equation (33) by F E £ yields 
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Since F T E d E by E c G and F T E d E by E c HE b E a Q are scalars, equation (34) becomes 
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The elements of the E and G matrixes are functions of the Laplace variable s. 
Equation (35) thus represents the transfer function of shock position to the perturbation 
(in total pressure for the example). By substituting s = jo>, the frequency response may 
be calculated by carrying out the matrix operations in complex algebra. This is done 
easily and quickly by a digital computer. Substitution of the shock position response into 
equation (33) as 

T 2 = HE b E a QP 0 + G^P 0 

*0 

gives the response for the elements of T 2 - Substituting these results into equation (27) 
gives the response for the elements of Tg and so on down the duct. Thus, the fre- 
quency response of total pressure, flow rate, and entropy are obtained at all stations. 

The response of static pressure at any station n is obtained using equation (23), which 
in matrix form is 
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COMPARISON OF ANALYSIS WITH METHOD OF CHARACTERISTICS 

Reference 11 gives transfer functions for shock position to perturbations in variables 
just ahead of the shock, and frequency response results that were compared with experi- 
mental data from reference 13. A method of characteristics solution was also included 



in reference 11 as a means for comparison since it is a very accurate analytical tech- 
nique (although computer time consuming). To evaluate the merit of the small perturba- 
tion analysis of this report as an analytical technique, it is compared in figure 4 with 
method-of-characteristics results from reference 11. The perturbed variable is Mach 
number just ahead of the normal shock. For the small perturbation analysis this is rep- 
resented as a flow perturbation with the aid of equation (59) from reference 9: 

AM = 5 + M 2 

M 5(M 2 - l) 

Assuming that area, total temperature, and pressure are held constant results in 



And in matrix form the perturbation equation (31) is 



Table I gives the numerical values for parameters used in the small perturbation 
analysis. No bleeds or bypasses were included in either analysis. 

The free-stream conditions were the same as the experimental test conditions (ref. 
13). The station Mach numbers are the same as those that were used in the method-of- 
characteristics solution and were supplied directly by the author. The A'/A param- 
eter calculation was based on the duct area distribution used in the method-of- 
characteristics analysis (ref. 11, fig. 28) and an assumed steady-state shock position 
41 centimeters aft of the cowl-lip station. 

The comparison of shock-position responses for the small perturbation and method 
of characteristics analyses (fig. 4) shows excellent agreement in both amplitude and 
phase. The results indicate that the small perturbation analysis is comparable in accu- 
racy to the method-of-characteristics analysis for this application. In comparing the 
results of either analysis with experimental results it should be remembered that neither 
technique includes viscous effects and that the problem is really three dimensional. 


^T + AA _ AW _ Aa r\ 
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Thus, the choice of coefficients, such as A'/A and Cg in the small perturbation 
analysis, will have considerable effect on the accuracy of the results. 

COMPARISON OF ANALYTICAL RESULTS WITH EXPERIMENTAL DATA 


Frequency response tests were conducted on a mixed-compression inlet in the Lewis 
10- by 10-foot supersonic wind tunnel. Perturbations both upstream and downstream of 
the normal shock were investigated and the results are reported in reference 13. As a 
means of verifying the analysis presented in this report, analytical results are compared 
with experimental data obtained from the upstream perturbations. 

A schematic diagram of the inlet and disturbance device setup is shown in figure 5. 
The disturbance was produced by a gust-generator plate located above and ahead of the 
inlet. The plate was oscillated sinusoidally at frequencies of 1 to 15 hertz to positive 
and negative angles of attack. This motion produces a perturbation in the flow-field 
Mach number and flow -angle ahead of the inlet. (A derivation for the disturbance will be 
given in the next section.) A frequency response of plate angle to command voltage is 
shown in figure 6. It shows that the amplitude decreases with increasing frequency and 
at 15 hertz is only 0.36 times its low-frequency amplitude. 

The inlet was sized for operation with a General Electric Company J85-13 turbojet 
engine and was alternately coupled to the engine or a long pipe. Provisions were made 
for choking the pipe airflow at either of two stations (fig. 5) and are designated as the 
short-pipe and long-pipe terminations. Design details of the inlet system and its steady- 
state performance characteristics are reported in reference 14, and a performance 
bleed investigation is reported in reference 15. 


Derivation of Gust Generator Airflow Perturbation 

Compression and expansion waves are generated by the leading edge of the plate 
when it is at positive and negative angles of attack (± cp), respectively. If only small and 
slow deflections in cp are considered, the first-order effect of <p on the static pres- 
sure ratio £ across the wave is (ref. 16, eq. (151)) 


I = 1 


+ 



1.4 Mj 



<P 


(38) 


where y=1.4 and <p is in radians. In terms of £ the Mach number Mg downstream 
of the wave is related to Mach number Mj upstream of the wave by (ref. 16, eq. (157)) 
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Linearization of equation (39) leads to (Mj = constant, 1=1, and Mj = M 2 when 
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From equation (38) A£ is again being constant) 


so that 


and since Mj = M 2 



(41) 


The first-order change in free -stream total pressure and entropy due to a compres- 
sion or expansion wave can be shown to be zero by the equations in reference 16. There- 
fore, 




= 0 


Neglecting the change in inlet capture area due to the change in free-stream Mach num- 
ber (calculated to be only about 0.4 percent per degree of plate angle for the experimen- 
tal inlet), equation (37) reduces to 


18 



or 
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Related to the plate angle perturbation the flow perturbation is (combining (41) and (42)) 
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And equation (31), relating the free-stream variables to the plate angle perturbation, is 
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Since the analysis is one-dimensional in nature, changes in flow angularity due to 
the plate angle of attack (limited to ±1°) were neglected. 

Analytical Computation Details 

A schematic of the inlet, to indicate the sectioning that was used for the analysis, is 
presented in figure 7. Comparisons of analytical and experimental frequency response 
results are shown in figures 8 to 10 with the inlet coupled to a turbojet engine, or a 
short-pipe or long-pipe termination. Except for the location of the choked exit station, 
the sectioning for each case was the same - all sections up to the exit being included. 
When coupled to the engine, the inlet was assumed to be choked at station 9. The short- 
pipe and long-pipe terminations were assumed to be choked at stations 10 and 12, 
respectively, which correspond to the actual choke points. 

The following information regarding sectioning is common to all three cases. The 
disturbance was assumed to occur at the intersection of the Mach wave generated by the 
plate at 0° angle-of-attack and the longitudinal axis of the inlet (see fig. 7). Choked 
bleed and bypass stations were located as shown in figure 7. The steady-state normal- 
shock position was assumed to be the same for all cases - the center of the throat 
static-pressure tap region shown in figure 4. Sections were included at the throat 


19 



exit-static pressure, static pressure before bypass cavity, and engine face static- 
pressure tap locations because it was desired to obtain frequency responses of those 
variables to the plate angle perturbation. The remaining sections were included to bet- 
ter approximate the variation in duct cross-sectional area. 

Tables II to IV contain numerical values for parameters used in the analysis, the 
results of which are presented in figures 8 to 10, respectively. The free-stream condi- 
tions correspond to the tunnel free-stream test conditions given in reference 13. These 
free-stream conditions plus the duct geometric cross-sectional area distribution were 
used as the basis for calculating the Mach number at each station, neglecting boundary- 
layer effects. The Mach number at the spike tip station was assumed to be the same as 
free stream. 

The value of the A'/A parameter was selected to make the analysis gain of shock 
position to plate angle of attack the same as the gain that was estimated from experimen- 
tal data. The values ranged from 0. 4534 to 1. 1453 per meter, depending on whether the 
inlet had stability bleed (reflected in the analysis by coefficient Cg). This range of 
values is approximately 3 to 7 times the value determined from the inlet geometric 
cross-sectional area variation at the shock operating point. Because of shock -boundary 
layer interaction effects, the shock gain is not uniquely dependent on the geometric 
value of A'/A. 

The steady- state airflow rate ratios W^W^, and and the 

stability bleed coefficient Cg were calculated with the aid of data found in references 
13 and 15, for the inlet-engine and inlet-long pipe cases. Unpublished NASA data were 
used to calculate these parameters for the inlet-short pipe case. 

DISCUSSION OF RESULTS 

Figure 8 shows comparisons of analysis and experimental (ref. 13) frequency 
response results with the inlet coupled to the J85-13 turbojet engine. 

In figures 8 to 10 two sets of analytical results are shown. The solid line results 
include all of the transportation times given in tables II to IV, respectively. For the 
dashed line results the transportation times in the seven supersonic sections were set 
equal to zero. This demonstrates the effect of neglecting dynamics in the supersonic 
portion of the inlet. 

The phase data are generally in good agreement for all signals. Even without the 
transportation times for the supersonic sections, the maximum difference in analysis 
phase results is about 15° (out of 100°) at 15 hertz. 

The amplitude data are presented as normalized amplitude ratios. Both experimen- 
tal and analytical amplitude ratio data for each response have been normalized by divid- 
ing by the low-frequency value of the amplitude ratio determined from experimental data. 
The amplitude data do not agree as well as the phase data, especially the static pressure 
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before bypass cavity and engine-face static-pressure responses. The analysis indicates 
that these two signals have low-frequency gains approximately 1.6 times larger than the 
experimentally measured gains. This discrepancy may be due in part to the fact that the 
engine-face station (compressor face) is not choked as was assumed for the analysis. 

The effect on analysis amplitude results is negligible when the supersonic section trans- 
portation times are eliminated. 

Comparisons of analytical and experimental (previously unpublished) results for the 
inlet-short pipe termination case are shown in figure 9. The plots for static pressure 
before the bypass cavity and for static pressure at the engine face are not given because 
these experimental responses were not taken. For this case both amplitude and phase 
data are in good agreement when the transportation times for the supersonic sections are 
included. Without the transportation times the difference in analysis phase is as much as 
25° (out of 100°) at 15 hertz. Shock position amplitude also shows more attenuation with- 
out the transportation times. 

Figure 10 shows a comparison of analytical and experimental (ref. 13) frequency 
response results for the inlet-long pipe combination. The plot for static pressure before 
the bypass cavity is not given because, again, this experimental response was not taken. 
As for the previous cases, the phase agreement is generally very good even with the 
transportation times for the supersonic sections set equal to zero. These times are less 
important in this case, because of the larger subsonic section transport times resulting 
from the longer duct. Agreement between the amplitude ratio data is good except for the 
throat-exit static -pressure signal. The analysis results indicate that the amplitude ratio 
is about 0.7 times the experimental low-frequency value. This error may be partly due 
to the fact that the analysis relates variables across the normal shock which ideally is a 
discontinuity. In nature, a normal shock may occur as a series of shocks or a shock 
train which could extend to the vicinity of the throat exit static pressure. 

Finally, a couple of observations will be made from the controls point of view re- 
garding the data in general. First, for each of the cases (figs. 8 to 10) the phase lag 
at a forward station is greater than at an aft station; or, the phase lag increases as 
signal distance from the perturbation location decreases. Second, comparing figures 
9 and 10 shows that the phase lag for a given signal is greater when the exit station is 
farther from the signal (long-pipe termination). These observations can also be made 
for the downstream disturbance case from data in reference 13. Thus, the inlet fre- 
quency response to a flow-field perturbation upstream of the normal shock is similar 
to the response resulting from a perturbation at the inlet exit. The normal shock must 
respond to the upstream perturbation before the subsonic duct pressures to satisfy the 
normal shock equations (7), (9), and (10). 
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CONCLUDING REMARKS 


An analysis was presented for determining the dynamic response of a supersonic 
inlet to perturbations ahead of the normal shock. The analysis, based primarily on the 
linearized, one -dimensional, distributed parameter wave equation, is applicable to analy- 
sis of mixed-compression inlets. Discontinuities and losses in the supersonic duct due 
to oblique shock waves are neglected. The equations are suitable for implementation on 
an analog computer. Also, closed-form expressions for the evaluation of frequency 
responses were obtained using matrix operations (no inversion required). These expres- 
sions are easily programmed on a digital computer and require little computer time to 
solve. 

The major difficulty in applying the analysis arises in choosing the correct value of 
the A’/A parameter. Because of shock-boundary layer interaction effects, the shock 
response is not uniquely dependent on the geometric value of A'/A. In this report an 
effective A'/A was determined empirically from the experimentally measured gain of 
shock position to the perturbation variable. The effective A'/A values varied from 3 
to 7 times the geometric value. Further work is required to determine a means for esti- 
mating an effective A'/A when experimental data are unavailable. This is important 
because A'/A is common to most inlet analyses. 

SUMMARY OF RESULTS 

Analysis frequency response results were compared with a method- of- characteris- 
tics solution over the frequency range of 5 to 40 hertz. The perturbed variable was Mach 
number just ahead of the normal shock. Agreement was excellent in both amplitude ratio 
and phase angle over the entire frequency range. The analysis has the advantages of 
being simpler to program and requiring much less computer time to solve. 

Frequency response comparisons were also made of the analysis with experimental 
inlet data obtained in a wind tunnel. The perturbation frequency range was 1 to 15 hertz. 
Data were obtained with the inlet coupled to a turbojet engine, a short pipe, or a long 
pipe (the pipes having choked exit airflows). In these cases the perturbed variable was 
Mach number just ahead of the inlet. 

Phase angle agreement was generally very good for both shock position and static 
pressures in the subsonic duct. When the transportation times for the supersonic duct 
were eliminated, the error in phase angle generally remained small except for the shock 
position responses with the inlet coupled to the engine or short pipe. 

The A'/A parameter (the rate of change of duct area with shock position divided 
by the duct area - evaluated at the shock operating point) was selected to assure good 
amplitude ratio agreement for the shock position responses. In general, the compari- 
sons showed good agreement between the shapes of the static-pressure amplitude ratio 
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response curves. However, in some cases a substantial shift in gain was observed. For 
the inlet-engine case the analysis predicted that static pressures near the subsonic dif- 
fuser exit had a low-frequency gain 1.6 times larger than the experimental values. The 
discrepancy may be due in part to the fact that the engine face is not choked as was 
assumed for the analysis. A static pressure near the normal shock was in error by as 
much as 27 percent. This error is believed to be due to the fact that the analysis relates 
variables across the normal shock which ideally is a discontinuity. In reality this shock 
may occur as a series of shocks or a shock train. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, September 9, 1974, 

501-24. 
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APPENDIX - SYMBOLS 



^1» C 4’ C 5 
C 2 ,C 3 

c Pn’ c Sn’ °Wn 
D 

E 

F 

G 

H 

i 

L 

M 

P 

P 

Q 

q 

R 

R c 

S 

s 

T 

t 

v 


2 

duct cross-sectional area, m 

2 

rate of change of duct cross-sectional area along duct, m /m 

speed of sound, m/sec 

coefficient matrix (3 by 3) 

stability bleed coefficient, dimensionless 

ratio of airflow rates at station 1 to most upstream station, dimension- 
less 

coefficients, dimensionless 
coefficients, sec 

coefficients at station n, dimensionless 
delay matrix (3 by 3) 
coefficient matrix (3 by 3) 
coefficient vector 

shock position transfer function vector 
coefficient matrix (3 by 3) 

V-T 

length, m 

Mach number, dimensionless 
o 

pressure, N/m 

Laplace transform of static pressure perturbation variable (see eq. (23)) 

coefficient vector 

elements of Q vector 

universal gas constant, J/mole-K 

radius of inlet cowl measured at cowl-lip station, m 

entropy, J/mole-K 

Laplace variable, sec -1 

state vector 

time, sec 

velocity, m/sec 
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W airflow rate, kg/sec 

X„ shock position, m 

9 

x space coordinate, m 

o', /3, y wave equation coefficients, dimensionless 

& perturbation quantity 

6 perturbation variable in Laplace domain 

e elements of E matrix 

d, a, r transportation times, sec 

| static pressure ratio across normal shock, dimensionless 

<p gust -generator plate angle of attack, rad or deg 

( t ) frequency, rad/sec 

Subscripts: 

a, b, c, d duct section identification for example problem 

av average value (one -half the sum of values at each end of duct section) 

B stability bleed 

by choked bypass or bleed 

E exit conditions 

1L 1L 

ij i in row, j in column location of matrix element 

n station or section number associated with variable 

s static condition 

T total or stagnation condition 

0, O’ stations upstream of normal shock for example problem 

1 station number just ahead of normal shock 

2 station number just aft of normal shock 

3, 4 stations downstream of normal shock for example problem 

00 most upstream station 

Superscripts: 

steady -state value of variable 

Laplace transform of nondimensional small perturbation variable 
T transpose of matrix 
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TABLE I.- NUMERICAL VALUES FOR PARAMETERS USED IN ANALYSIS COMPARISON 
WITH METHOD-OF-CHARACTERISTICS SOLUTION 3 

[Free-stream test conditions: total pressure, 10.05 N/cm 2 ; total temperature, 342 K; 
Mach number, 2.497; area ratio parameter, A'/A, 0.961 m"*!] 


Station 

Distance 

from 

cowl-lip, 

cm 

Mach 

number 

Section 

between 

stations 

Transportation time, msec 

(T 

T 

e 

1 

41 

1.3 


— 

— 

— 

2 

41 

.786 

lto 2 

— 


— 

3 

47 

.684 

2 to 3 

0.098 

0.642 

0.232 

4 

54 

.651 

3 to 4 

.118 

.592 

.295 

5 

60 

.592 

4 to 5 

.104 

.444 

.270 

6 

70 

.479 

5 to 6 

.181 

.597 

.519 

7 

80 

.427 

6 to 7 

.189 

.503 

.608 

8 

90 

.391 

7 to 8 

.195 

.464 i 

.670 

9 

100 

.363 

8 to 9 

.199 

.439 

.725 

10 

110 

.369 

9 to 10 

.200 

.431 

.747 

11 

120 

.381 

10 to 11 

.199 

.437 

.730 

12 

130 

.369 

11 to 12 

.199 

.437 

.730 

13 

140 

.358 

12 to 13 

.200 

.429 

.752 


a See fig. 4. 
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TABLE H. - NUMERICAL VALUES FOR PARAMETERS USED IN INLET WITH ENGINE ANALYSIS 3 

[Free-stream test conditions: total pressure, 9.98 N/cm^; total temperature, 342. OK; 

Mach number, 2.498. Stability bleed coefficient^ 0; a_rea ratio parameter, A'/A, 

0.4534 m"*. Steady-state airflow rate ratios: /w^ at station 5’, 0.015; 

W by /w M at station 3', 0.025; W byy /w B at station 6, 0.116; Wj/w m , 0.958; 
at station 9, 0.842] 


Station 

Mach 

number 

Section 

between 

stations 

Transportation time, msec 

a 

T 

e 

S' 

2.498 


— 


— 

V 

2.498 

8 f to V 

0.793 

-1.851 

1.110 

6 T 

2.252 

V to 6 T 

.554 

-1.363 

.788 

5 T 

1.513 

6 ? to 5 T 

.359 

-1.201 

.553 

4 f 

1.539 

5 f to 4 T 

0 

0 

0 

3 f 

1.383 

4 f to 3 f 

.046 

- .248 

.078 

2' 

1.434 

3’ to V 

0 

0 

0 

1 

1.390 

V to 1 

.116 

- .677 

.198 

2 

.744 

1 to 2 

0 

0 

0 

3 

.722 

2 to 3 

.094 

.605 

.221 

4 

.640 

3 to 4 

.139 

.729 

.342 

5 

.388 

4 to 5 

.676 

2.096 

1.998 

6 

.377 

5 to 6 

.305 

.682 

1.102 

7 

.324 

6 to 7 

0 

0 

0 

8 

.335 

7 to 8 

.456 

.905 

1.842 

9 

.322 

8 to 9 

.106 

.210 

.430 


^ee fig. 8. 
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TABLE HI. - NUMERICAL VALUES FOR PARAMETERS USED IN INLET WITH SHORT -PIPE -ANALYSIS 4 

[Free-stream test conditions: total pressure, 8.94 N/cm 2 ; total temperature 313. 4 K; Mach 
number, 2.498. Stability bleed coefficient, 0.037; area ratio parameter, A’/A, 1.1453 m* 1 . 
Steady-state airflow rate ratios: W^y/Woo at station 5', 0.015j_ ^by/^°° a * station 3', 

0.025; at station 6, 0.135; 0.958; W E/ /W w at station 10, 0.806] 


Station 

Mach 

number 

Section 

between 

stations 

Transportation time, msec 


a 

r 

e 

8 f 

7* 

2.498 

2.498 

8 T to 

V 

0.828 

-1.934 

1.160 

6* 

2.252 

7' to 

6* 

.579 

-1.424 

.823 

5* 

1.514 

6 T to 

5* 

.375 

-1.255 

.578 

4* 

1.539 

5’ to 

4' 

0 

0 

0 

3 1 

1.383 

4’ to 

3' 

.048 

- .259 

.081 

V 

1.435 

3’ to 

V 

0 

0 

0 

1 

1.407 

2’ to 

1 

.121 

- .694 

.205 

2 

.737 

1 to 

2 

0 

0 

0 

3 

.698 

2 to 

3 

.098 

.597 

.235 

4 

.621 

3 to 

4 

.146 

.713 

.368 

5 

.377 

4 to 

5 

.712 

2.123 

2.145 

6 

.365 

5 to 

6 

.321 

.699 

1.185 

7 

.305 

6 to 

7 

0 

0 

0 

8 

.313 

7 to 

8 

.483 

.917 

2.045 

9 

.301 

8 to 

9 

.113 

.213 

.479 

10 

.345 

9 to 10 

.362 

.708 

1.482 


a See fig. 9. 
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TABLE IV. - NUMERICAL VALUES FOR PARAMETERS USED IN INLET WITH LONG-PIPE -ANALYSIS 51 

(Free-stream test conditions: total pressure, 9.02 N/cm^; total temperature, 316. OK; 

Mach number, 2.498. Stability bleed co efficient^ 0; area ratio parameter, A '/A, 

JX4716m“*. Steady-state airflow rate ratios: ^by/^°° at station 5’, 0.015; 

W by /_W M at station 3’, 0.025; W by/ /w M at station 6, 0.098; W^W^, 0.958; 

^E/"^oo at station 12, 0.860] 


Station 

Mach 

number 

Section 

between 

Transportation time, 

msec 


stations 

a 

T 

e 

8 f 

2.498 





V 

2.498 

8’ to V 

0.825 

-1.926 

1.155 

6 1 

2.252 

V to 6' 

.577 

-1.418 

.820 

5' 

1.513 

6 T to 5 f 

.374 

-1.250 

.575 

4* 

1.538 

5' to 4' 

0 

0 

0 

3’ 

1.382 

4 T to 3 f 

.048 

- .258 

.081 

2* 

1.434 

3 f to V 

0 

0 

0 

1 

1.390 

V to 1 

.120 

- .705 

.206 

2 

.744 

1 to 2 

0 

0 

0 

3 

.718 

2 to 3 

.097 

.626 

.230 

4 

.634 

3 to 4 

.144 

.747 

.358 

5 

.379 

4 to 5 

.707 

2.147 

2.107 

6 

.366 

5 to 6 

.319 

.698 

1.176 

7 

.323 

6 to 7 

, o 

0 

0 

8 

.330 

7 to 8 

.476 

.937 

1.934 

9 

.316 

8 to 9 

.111 

.217 

.455 

10 

.362 

9 to 10 

.357 

.723 

1.408 

11 

.330 

10 to 11 

.604 

1.243 

2.348 

12 

.342 

11 to 12 

4.161 

8.375 

16.536 


a See fig. 10. 


Y (1) Conditions just upstream 
\of shock 
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Main duct 
airflow 


Bleed or bypass 
airflow 

1 1 


[- — Zero length assumed 


I 


Station: n n+1 

Figure 2. - Bleed or bypass airflow in inlet duct. 


Station q o' 1 2 3 4 E 



Figure 3. - Schematic of inlet for example problem. 




Figure 4. - Comparison of small perturbation analysis 
with method-of-characteristics solution. Response 
of shock position to perturbation in Mach number 
just upstream of shock. Normalizing amplitude 
ratio, 2.067. (Numerical values of parameters used 
irxanaJysis given in table I. > 


Throat-exit static 
pressure 
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Phase angle, deg Normalized 

amplitude ratio 



Frequency, Hz 


(a) Signal, shock position; (b) Signal, throat-exit static (c) Signal, static pressure (d) Signal, engine face 
normalizing amplitude pressure; normalizing before bypass cavity; nor- static pressure; normal- 

ratio, (AX s /R c )/A<p, amplitude ratio, malizing amplitude ratio, izingamplitude ratio, 

0.3381 deg-1. (AP/P)/A<p, 0.0464 deg -1 . (AP/P)/A<P, 0.0245 deg -1 . (AP/P)/A<P, 0.0239 deg -1 . 

Figure 8. - Comparison of analysis with experimental inlet-engine response to gust-generator plate. (Numerical values 
for parameters used in analysis given in table II. ) 



Frequency, Hz 

(a) Signal, shock position; normal- (b) Signal, throat exit static pressure; 

izing amplitude ratio, normalizing amplitude ratio, 

(AX s /R c )/A<0, 0.1039 deg" 1 . (AP/P)/Ag>, 0.0375 deg" 1 . 

Figure 9. - Comparison of analysis with experimental inlet-short-pipe response to 
gust-generator plate. (Numerical values for parameters used in analysis given 
in table III.) 
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